Observation of Z) sl (2536) + — ► D + n K + and angular 
decomposition of D al (2S36) + — > D*+K% 

V. Balagura, 14 I. Adachi, 9 H. Aihara, 44 K. Arinstein, 1 V. Aulchenko, 1 T. Aushev, 19 ' 14 

A. M. Bakich, 40 E. Barberio, 22 A. Bay, 19 K. Belous, 13 V. Bhardwaj, 34 U. Bitenc, 15 
A. Bondar, 1 A. Bozek, 28 M. Bracko, 21 ' 15 J. Brodzicka, 9 T. E. Browder, 8 M.-C. Chang, 4 
Y. Chao, 27 A. Chen, 25 W. T. Chen, 25 B. G. Cheon, 7 R. Chistov, 14 I.-S. Cho, 49 Y. Choi, 39 
J. Dalseno, 22 M. Danilov, 14 A. Drutskoy, 3 S. Eidelman, 1 D. Epifanov, 1 N. Gabyshev, 1 

A. Garmash, 35 B. Golob, 20 ' 15 H. Ha, 17 J. Haba, 9 T. Hara, 33 N. C. Hastings, 44 
K. Hayasaka, 23 H. Hayashii, 24 M. Hazumi, 9 D. Heffernan, 33 Y. Hoshi, 42 W.-S. Hou, 27 
H. J. Hyun, 18 K. Inami, 23 A. Ishikawa, 36 H. Ishino, 45 R. Itoh, 9 M. Iwasaki, 44 Y. Iwasaki, 9 
D. H. Kah, 18 H. Kaji, 23 N. Katayama, 9 H. Kawai, 2 T. Kawasaki, 30 H. Kichimi, 9 
H. J. Kim, 18 H. 0. Kim, 39 Y. J. Kim, 6 K. Kinoshita, 3 S. Korpar, 21 * 15 P. Krizan, 20 - 15 
P. Krokovny, 9 C. C. Kuo, 25 A. Kuzmin, 1 Y.-J. Kwon, 49 J. S. Lange, 5 M. J. Lee, 38 
S. E. Lee, 38 T. Lesiak, 28 A. Limosani, 22 S.-W. Lin, 27 Y. Liu, 6 D. Liventsev, 14 F. Mandl, 12 
S. McOnie, 40 T. Medvedeva, 14 W. Mitaroff, 12 H. Miyake, 33 H. Miyata, 30 Y. Miyazaki, 23 
R. Mizuk, 14 D. Mohapatra, 48 G. R. Moloney, 22 Y. Nagasaka, 10 E. Nakano, 32 M. Nakao, 9 
H. Nakazawa, 25 Z. Natkaniec, 28 S. Nishida, 9 0. Nitoh, 47 S. Ogawa, 41 T. Ohshima, 23 
S. Okuno, 16 S. L. Olsen, 8 - 11 W. Ostrowicz, 28 H. Ozaki, 9 P. Pakhlov, 14 G. Pakhlova, 14 
H. Palka, 28 C. W. Park, 39 L. S. Peak, 40 R. Pestotnik, 15 L. E. Piilonen, 48 Y. Sakai, 9 
O. Schneider, 19 K. Senyo, 23 M. Shapkin, 13 C. P. Shen, 11 H. Shibuya, 41 J.-G. Shiu, 27 

A. Somov, 3 S. Stanic, 31 M. Staric, 15 T. Sumiyoshi, 46 K. Tamai, 9 M. Tanaka, 9 
G. N. Taylor, 22 Y. Teramoto, 32 I. Tikhomirov, 14 S. Uehara, 9 K. Ueno, 27 T. Uglov, 14 
Y. Unno, 7 S. Uno, 9 P. Urquijo, 22 Y. Usov, 1 G. Varner, 8 K. Vervink, 19 S. Villa, 19 
A. Vinokurova, 1 C. C. Wang, 27 C. H. Wang, 26 M.-Z. Wang, 27 P. Wang, 11 X. L. Wang, 11 

Y. Watanabe, 16 E. Won, 17 B. D. Yabsley, 40 A. Yamaguchi, 43 Y. Yamashita, 29 
M. Yamauchi, 9 Z. P. Zhang, 37 V. Zhilich, 1 V. Zhulanov, 1 A. Zupanc, 15 and O. Zyukova 1 

(The Belle Collaboration) 

1 Budker Institute of Nuclear Physics, Novosibirsk 

2 Chiba University, Chiba 
3 University of Cincinnati, Cincinnati, Ohio 45221 
^Department of Physics, Fu Jen Catholic University, Taipei 
5 Justus-Liebig-Universitat Giefien, Giefien 
6 The Graduate University for Advanced Studies, Hayama 
7 Hanyang University, Seoul 
8 University of Hawaii, Honolulu, Hawaii 96822 
9 High Energy Accelerator Research Organization (KEK), Tsukuba 
10 Hiroshima Institute of Technology, Hiroshima 
11 Institute of High Energy Physics, 
Chinese Academy of Sciences, Beijing 
12 Institute of High Energy Physics, Vienna 
13 Institute of High Energy Physics, Protvino 
1Jt Institute for Theoretical and Experimental Physics, Moscow 
15 J. Stefan Institute, Ljubljana 



Typeset by REVTeX 



1 



16 Kanagawa University, Yokohama 
17 Korea University, Seoul 
18 Kyungpook National University, Taegu 
19 Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne 
20 University of Ljubljana, Ljubljana 
21 University of Maribor, Maribor 
22 University of Melbourne, School of Physics, Victoria 3010 
23 Nagoya University, Nagoya 
24 Nara Women's University, Nara 
25 National Central University, Chung-li 
26 National United University, Miao Li 
21 Department of Physics, National Taiwan University, Taipei 
28 H. Niewodniczanski Institute of Nuclear Physics, Krakow 
29 Nippon Dental University, Niigata 
30 Niigata University, Niigata 
31 University of Nova Gorica, Nova Gorica 
32 Osaka City University, Osaka 

33 Osaka University, Osaka 
3Jl Panjab University, Chandigarh 
35 Princeton University, Princeton, New Jersey 08544 
36 Saga University, Saga 
37 University of Science and Technology of China, Hefei 
38 Seoul National University, Seoul 
39 Sungkyunkwan University, Suwon 
4° University of Sydney, Sydney, New South Wales 
41 Toho University, Funabashi 
42 Tohoku Gakuin University, Tagajo 
^ 3 Tohoku University, Sendai 
^Department of Physics, University of Tokyo, Tokyo 

45 Tokyo Institute of Technology, Tokyo 

46 Tokyo Metropolitan University, Tokyo 

^ Tokyo University of Agriculture and Technology, Tokyo 
48 Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 

^Yonsei University, Seoul 



2 



Abstract 

Using 462 fb" 1 of e + e~ annihilation data recorded by the Belle detector, we report the first ob- 
servation of the decay D s i(2536) + — > D + n~K + . The ratio of branching fractions 13< ^^+ >D D *+^o) " 
is measured to be (3.27 ± 0.18 ± 0.37)%. We also study the angular distributions in the 
-D s i(2536) + — ► D* + Kg decay and measure the ratio of D- and S-wave amplitudes. The S-wave 
dominates, with a partial width of rs/r tota j = 0.72 ± 0.05 ± 0.01. 

PACS numbers: ll.80.Et, 13.25.Pt, 13.30.Eg, 13.66.Bc, 13.87.Fh, 13.88. +e, 14.40.Lb 
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INTRODUCTION 



Two states, -D* (2317) + and _D s i(2460) + , have been discovered recently both in continuum 
e + e~ annihilation near ^/s = 10.6 GeV/c 2 and in B meson decays [H, 0, [1] . Their spin-parities 
are, respectively, J p = + and 1 + (4j, and they are presumed to be P-wave excited cs states 
with j = L + S s = 1/2. Here, L = 1 is the orbital angular momentum and S s is the spin 
of the light antiquark. However, their masses are unexpectedly low |5|]. This has renewed 
interest in measurements of P-wave excited charm mesons. 

We report the first observation of the decay _D s i(2536) + — > D + 7r~K + . (The inclusion of 
charge-conjugate modes is implied throughout this paper.) The D + n~ pair in the final state 
is the only Dir combination that cannot come from a D* resonance: D*° mesons can only 
be produced virtually here since M D *o < M D + + M % -. The new Z) s i(2536) + — > D + n~K + 
mode reported here is only the second observed three-body decay of the D s i(2536) + , after 

D+7T+7T- Q. 

In addition, we have performed an angular analysis of the _D s i(2536) + — > D* + Kg mode. 
Heavy Quark Effective Theory (HQET) predicts that for an infinitely heavy c-quark this 
decay of a J p = 1 + , j = 3/2 state should proceed via a pure D-wave p. The corresponding 
decay of its partner, the _D s i(2460) + , which is believed to be a 1 + , j = 1/2 state is energetically 
forbidden, but if it were allowed it would proceed via a pure S-wave. Since heavy quark 
symmetry is not exact, the two 1 + states can mix with each other, 



|A,i(2460)+) = cos6 l 1 / 2 ^ + sin6 | 3 / 2 £i), 

|D s i(2536)+) = -sin^ \ 1 / 2 E 1 ) + cose \V 2 E^ ' ' 



where \ l ^ 2 Eij and \ z / 2 E\j denote the states with j = 1/2 and j = 3/2, respectively. Note that 
the coupling via common decay channels can give a contribution to the mixing that might not 
be well represented by an orthogonal rotation [7| . We neglect this possibility in the expression 
above. If 9 ^ 0, an S-wave component can appear in the decay Z) s i(2536) + — > D*K. 
Moreover, even if 6 is small, the S-wave component can give a sizeable contribution to the 
width because the D-wave contribution is strongly suppressed by the small energy release 
in the D,i(2536)+ -> D*K decay. 

The first attempt to decompose S- and D-waves in the analogous decays of the non- 
strange mesons _D!(2420)° — >• D* + n~ and _D!(2420) + — > D*°n + was reported more than 
ten years ago by CLEO @, Q; currently, no results on the _D s i(2536) + exist. Moreover, 
CLEO's method did not allow the measurement of the ratio of partial widths: it only 
determined the relation between this ratio and the relative phase between the S- and D- 
wave amplitudes. Some information on e is obtained from the ratio of electromagnetic decay 
rates £> s i(2460) + — ► £>+7, -D* + 7, since only the 1 P± state in D sl (2A60) + undergoes an El 
transition to _D+ and only the 3 P 1 state to D* + [Io|. The bases \ j Ei) and \ 2S+1 P^ are 
related by the rotation angle 9 , where tan^o = — \/2. The angle between the bases \Df) 
and | 25+1 Px) is 6 + 9 . The Belle Collaboration studied D sl {2460) + -> D+7, D* + ^ decays 
using D sl (2460) + from both B decays [2| and from e + e~ annihilation 0, and determined 



the ratio of decay rates to be 0.4 ± 0.3 and 0.28 ± 0.17, respectively. According to Ref. [10 
the average ratio of 0.31 ± 0.14 gives the constraint tan 2 (6 l + 9 Q ) = 0.8 ± 0.4. Detailed 
knowledge of the mixing is important to test different theoretical models 0, GU 11 1 . to fix 
their parameters and to understand better the nature of D s j mesons. 
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Finally, using the Z) s i(2536) + — > D* + K s mode, we have measured the spin alignment of 
high momentum D s i(2536) + mesons produced in e + e~ annihilation. Production of excited 



mesons in the HQET framework is described in Ref. [12|. The fragmentation process is 
assumed to be so rapid that the color magnetic forces do not have time to act and thus the 
spin of the light antiquark in the produced meson is uncorrelated with that of the heavy 
quark. One consequence of this is that D* mesons with j = \ are produced unpolarized. This 
was confirmed with good accuracy by CLEO [l3[ in e + e~ — > cc events at yfs = 10.5 GeV and 



was also checked by other experiments [14|. Another prediction is that D* and D mesons 
are produced according to the number of available helicity states in a 3:1 ratio. However, 
experimental data from several different production mechanisms (e + e~, hadropro duct ion, 
photoproduction, etc.) give an average probability for an S-wave meson to be produced in 
a vector state of 0.594 ± 0.010 [15], which is much smaller than the expected value of 0.75. 

There are no similar measurements for the P-wave states. Contrary to the (D, D*) 
case, HQET predicts that the members of j = 3/2 doublet can be produced aligned. The 
probabilities for the light degree of freedom to have helicity —3/2, —1/2, 1/2, 3/2 are 
expressed via one parameter w^/ 2 as §?£>3/2, |(1 — ^3/2), |(1 — ^3/2), ^ w 3/2, respectively. 
By adding the c-quark spin and resolving the cs system into 1 + and 2 + states, one can 
calculate their alignment. For D sl (2536) + the probability of zero helicity is p o — ^3/2)- 
A calculation based on perturbative QCD and a nonrelativistic quark model gives w%/ 2 = 
29/114 m 0.254 [l6[ and poo ~ 0.497. This calculation also predicts the dependence of w 3 / 2 on 
the longitudinal momentum fraction and on the transverse momentum of the meson relative 
to the heavy quark jet. The ARGUS analysis of the angular distributions in 1)2(2460) — > Dn 
decay [l7[ gives an upper limit w 3 / 2 < 0.24 at 90% CL [12]. Once w 3 / 2 is measured, one 
can make definite predictions for the angular distributions of the remaining j = 3/2 meson 
decays and check the validity of HQET. 



SELECTION CRITERIA 



This study is based on a data sample of 462 fb 1 collected near the T(4S) resonance with 



the Belle detector at the KEKB asymmetric-energy e + e~ (3.5 on 8 GeV) collider [18j. The 
Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon vertex 
detector, a 50-layer central drift chamber, an array of aerogel threshold Cherenkov counters, 
a barrel-like arrangement of time-of-flight scintillation counters, and an electromagnetic 
calorimeter comprised of CsI(Tl) crystals located inside a superconducting solenoid coil that 
provides a 1.5 T magnetic field. An iron flux-return located outside the coil is instrumented 
to detect K\ mesons and to identify muons. The detector is described in detail elsewhere [l9| . 
Two inner detector configurations were used. A 2.0 cm beampipe and a 3-layer silicon vertex 
detector were used for the first sample of 155 fb^ 1 , while a 1.5 cm beampipe, a 4-layer silicon 
detector and a small-cell inner drift chamber were used to record the remaining 307 ftr 1 [20]. 

K and ^ candidates are required to originate from the vicinity of the event-dependent 
interaction point. To identify kaons, we combine the ionization energy loss (dE/dx) from 
the central drift chamber, time of flight and Cherenkov light yield information for each 
track to form kaon and pion likelihoods Ck and C n , respectively, [2l| and then impose 
the requirement £k/(£k + £-k) > 0.1. This requirement has 98% (97%) efficiency for a 
kaon from D sl (2536) + (kaon from D) and a 12% (17%) misidentification probability for a 
pion with the same momentum. All unused tracks, whether identified as a kaon or not, 
are treated as pion candidates in what follows. K s candidates are reconstructed via the 
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7r + 7r~ decay channel, with a mass within ±30MeV/c 2 of the nominal K s mass (among 
other quality requirements). D° and D + mesons are reconstructed using K~tt + , K s it + 7t~, 
K~7t + 7t + 7t~ and K s ir + , K~7t + tt + decay modes, respectively. All combinations with masses 
within ±20 MeV/c 2 of the nominal D mass are selected (99% efficiency); a mass and vertex 
constrained fit is then applied. 

D* + mesons are reconstructed using the D tt + mode. The slow 7r + momentum resolution 
is degraded by multiple scattering, but is improved by a track refit procedure in which the 
7T+ origin point is constrained by the intersection of the D° momentum and the known e + e~ 
interaction region. The D°tt + mass is required to be within ±1.5 MeV/c 2 of the D* + nominal 
value, which corresponds to 98% efficiency A D* + mass constraint is not imposed. Instead, 
we characterize the _D s i(2536) + candidate using the mass difference M D o n + K o —M D o n +, where 
the error in the D* + momentum nearly cancels out. For the -D s i(2536) + — > D + ix~K + decay 
mode, the track refit procedure described above is applied to the pion and kaon momenta, 
and the Z) s i(2536) + is characterized by the D + 7r~K + mass. 

It is known that the momentum spectrum of the excited charm resonances from continuum 
e + e~ annihilation is hard. In addition, due to the strong magnetic field in the Belle detector, 
the reconstruction efficiency for slow 7r =1= and K + mesons rises with D sl (2536) + momentum. 
Therefore, we require xp > 0.8 for the scaled momentum xp, defined as the ratio p* /p* max - 
Here, p* is the momentum of the -D s i(2536) + candidate in the e + e~ center-of-mass frame, 

while Pmax = \l -^beam — M 2 is the momentum in this frame for a candidate carrying all the 
beam energy. This selection also removes D sl (2536) + mesons produced in the decays of B 
mesons. 

In Monte Carlo (MC) simulation, D sl (2536) + mesons from e + e~ annihilation, particle 
decays and the detailed detector response are simulated using the PYTHIA, EvtGen and 
GEANT packages [22I], respectively. The D° and D + decay modes used in reconstruction 
are generated with their resonant substructures taken from the Particle Data Group (PDG) 
compilation [4j but neglecting any interference effects. The _D s i(2536) + momentum spectrum 
as measured with the _D s i(2536) + — > D* + K s decay mode is used for MC generation. As 
shown below, no clear resonant substructure is visible in the decay _D s i(2536) + — ► D + n~K + . 
Therefore, this mode is simulated as a three-body phase space decay. _D s i(2536) + — > D* + K s 
decays are generated according to our measured _D s i(2536) + polarization and D/S-wave 
interference. 

L» s i(2536)+ -» D + ir-K + DECAY 

The mass M D + n - K + (upper plot) and the mass difference (M^o^+^o — M D o w +) + M™+ 
(lower plot) for all accepted combinations are shown in Fig. [lj The PDG superscript denotes 
the nominal mass value from Ref. [4[]. A clear peak for the new decay channel _D s i(2536) + — ► 
D + 7i~K + is visible. The mass spectrum of the wrong sign combinations D + tc + K~ shown 
by the hatched histogram has no enhancement in the _D s i(2536) + region. 

To obtain the number of -D s i(2536) + decays, each of the distributions in Fig. [T] is fit 
to the sum of two Gaussians with a common mean (but not necessarily common between 
the two decay modes). To ensure that the second Gaussian is always wider than the first 

one, its width is chosen to be of the form a 2 = y 1 <j\ ± Act 2 . The position of the peak, 
01, Act, the fraction of events in the first Gaussian and the total number of events in two 
Gaussians are allowed to vary in the fit. The background for the three-body D + n~K + 
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2500 



2510 2520 2530 2540 2550 2560 
M(D 7i + Ks) - M(DV) + M(Dpdg) MeV/c 2 



FIG. 1: D s i(2536)+ mass spectra for D + tt-R + (top) and D* + K s (bottom) decay modes. The 
hatched histogram in the top plot shows the corresponding spectrum of wrong sign D + n + K~ 
combinations. The fit is described in the text. The fit results are listed in Table HI 



Kg) mode is parameterized by a second (first) order polynomial multiplied 

2 



(two-body D 

by the threshold function (M - M^ 7r _ it+ ) 2 (^M - M^ +k0 ), where Mf r is the sum of the 

nominal masses of final state particles / Q. TableU contains the fit results together with 
the parameters of the Gaussians obtained from MC simulation. There is a small fraction 
of events that contribute two entries to the D s i(2536) + signal region in the mass plot. The 
excess of such events in comparison with the same number averaged over the left and the 
right sideband is 35-18.5 = 16.5 and 203-15.5 = 187.5 for the D + ti~K + and D*+K s modes, 
respectively. These values are estimates of the contribution of double counted signal events. 
In the branching ratio calculation they are subtracted from the yields given in Table [H The 
signal and the sidebands are denned as \ AM D +J < 5 MeV/c 2 , 10 MeV/c 2 < \AM D + \ < 20 
MeV/c 2 , respectively, where AM D + is measured relative to the peak position obtained from 
the fit. 

The ratio of branching fractions is found to be 



£(£> s i(2536)+ -> D+TT-K+) 
£(£>,! (2536)+ ->• D*+K°) 



(3.27 ±0.18 ±0.37)%, 



(2) 



where the first error is statistical and the second is systematic. This ratio takes into account 
the partial branching fractions for mesons to the selected final states and of neutral 
kaons to 7r + 7r~. 

The systematic error receives contribution from the sources listed in Table [TTJ A possible 
difference between the data and MC simulation in evaluation of the tracking efficiency is 
estimated using partially reconstructed D* + decays. The reconstruction efficiency errors 
for slow K + , 7r ± and K s are added linearly. The uncertainty in the kaon particle identifi- 
cation is estimated using D* + decays. Uncertainty in the ratio of D + and D° efficiencies 
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TABLE I: Fit results for the -D s i(2536) + spectra in Fig. [T] and for the corresponding MC sim- 
ulation spectra: number of events in the two Gaussians for data (or the efficiency for MC sim- 
ulation), fraction of events in the narrow Gaussian, width of the narrow Gaussian, additional 
width contribution for the wide Gaussian, and the mass difference with respect to Mp® G = 
(2535.35 ± 0.34 ± 0.5) MeV/c 2 . 





Yield 


Narrow Gaussian 


0"l 


Ao- 






(Efficiency) 


Fraction 


(MeV/c 2 ) 


(MeV/c 2 ) 


(MeV/c 2 ) 


D + n~K + , data 


1281 ±66 


0.59 ±0.06 


0.76 ±0.06 


2.4 ±0.4 


-0.57 ±0.04 


D + it-K + , MC 


0.2699 ± 0.0017 


0.463 ± 0.014 


0.94 ±0.02 


2.55 ± 0.04 


-0.031 ±0.010 


D* + K s , data 


5673 ± 81 


0.63 ±0.03 


1.01 ±0.03 


2.54 ±0.13 


-0.43 ±0.02 


D* + K s , MC 


0.1273 ± 0.0004 


0.629 ± 0.006 


0.946 ± 0.008 


2.56 ±0.03 


-0.034 ±0.005 



is conservatively estimated by a comparison of different decay modes used in the recon- 
struction. One of the largest contributions to the systematic uncertainty arises due to the 
model of the background. It is evaluated by fitting the wrong sign D + tt + K~ subtracted 
spectrum, which contains almost no background. For the D + ir~K + mode, the efficiency 
is almost independent of the D + ir~ , K + tt~ masses and the angular distribution of decay 
products. Therefore, the possible difference between the simplified phase space MC model 
and the real D s i(2536) + — ► D + n~K + decay results in a small uncertainty in the efficiency 
determination. It is estimated by comparing the yields of events using either an average or 
differential efficiency in the D sl (2536) + decay angles and the D + n~ and K + w~ masses. The 
total systematic error is found to be 11.2% (Table [IB. 



TABLE II: Systematic uncertainties for B(D s i(2536)+ -» D + tT K + )/B(D s i(2536)+ -> D* + K°). 

Source Uncertainty, % 

Reconstruction efficiencies of slow tt^, K + and K s from D* + and D s i(2536) + 7.5 

Particle identification of slow K + from Z? s i(2536) + 1.2 

Ratio of D + and D° efficiencies 2.7 

Background model in M (D + ir~ K + ) spectrum 6.5 
Efficiency dependence on D + ir~ , K + ir~ masses and 

angular distribution of decay products in D + ir~ K + decay 1.2 

Branching ratio of intermediate resonances [4] 4.1 

Total 11.2 



To cross-check the results, the D + mass spectrum is plotted in Fig. [2] for the D s i(2536) + 
signal and sidebands. The latter is normalized to the area of the signal interval. The sideband 
subtracted plot shown in the bottom of Fig. [2] is fit to a double Gaussian as above and a 
constant background. The resulting yield 1249 ± 66 is consistent with the yield 1262 ± 65 
obtained from the fit of the D s i(2536) + mass spectrum. The constant background level is 
found to be —0.9 ± 0.8, which is consistent with zero. The enhancement in the D + mass 
region observed in the _D s i(2536) + sidebands (top plot of Fig. [2]) is due to combinations of 
a real D + with a random ir~K + pair in the event. 

The D + 7i~ and K + tc~ mass distributions for the _D sl (2536) + — ► D + -k~K + decay are 
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FIG. 2: D + mass spectrum for the D s i(2536) + signal band (\AM D + % - K +\ < 5 MeV/c 2 , open 
histogram in the top plot) and the sidebands (10 < | AM D + n - K +\ < 20 MeV/c 2 , normalized to the 
signal interval, hatched histogram). AM D + 1T - K + is measured relative to the peak position in the 
top plot of Fig. [TJ The bottom plot shows the sideband subtracted distribution. The solid curve 
shows the results of the fit described in the text. 



shown in Fig. [3j The -D s i(2536) + signal yield is obtained from fits to the D + n~ K + mass 
distribution in bins of D + n~ and K + tt~ mass. All Gaussian parameters except the total 
number of events are fixed in the fit to the values listed in Table [H The position of the 
threshold used for the background description depends on the chosen bin. The areas under 
the histograms have been normalized to unity. The spectra are not efficiency corrected. The 
dashed histograms show the corresponding MC spectra for D sl (2536) + — > D + -k~K + decays 
simulated according to a phase space distribution. From this plot, the data points do not 
appear to be entirely consistent with a phase space distribution, but neither do they exhibit 
any clear dominant resonant substructure. 



ANGULAR ANALYSIS OF D s i(2536)+ D*+K s DECAY 

The _D s i(2536) + — > D* + K s decay kinematics can be described by three angles a, (3 and 7 
defined as shown in Fig. |H The angles a and /3 are measured in the rest frame: a is the 
angle between the boost direction of the e + e~ center-of-mass and the Kg momentum, while 
(3 is the angle between the plane formed by these two vectors and the decay plane. The 
third angle 7 is defined in the D* + rest frame between 7r + and Kg. 

/} s i(2536) + polarization can be described in terms of its helicity density matrix p mim2 . 
The contribution of the element p mi m 2 to the decay amplitude is proportional to e ~ l ^ mi ~ m2 \ 
where <fi is the azimuthal rotation angle around the e + e~ boost direction in the rest 
frame. After integration over the contribution of off-diagonal elements vanishes. Due to 
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FIG. 3: Normalized mass spectra of D + ir~ (top) and K^tt~~ (bottom) pairs from D s i(2536) + — 
D + ir~ K + decay obtained from fits to the D + tt~ K + mass distributions in different D + ir~ or K + ir~ 
mass bins. The dashed histograms show the corresponding MC distributions for D sl (2536) + — 
D + ir~K + decays simulated according to a phase space distribution. 



Boost 




FIG. 4: Definitions of the angles a, [3 and 7. The first two are measured in the D s i(2536) + rest 
frame, the third in the D* + frame. "Boost" refers to the direction of the e + e~ center-of-mass in 
the -D s i(2536) + rest frame. 



parity conservation, the three diagonal elements can be expressed in terms of the longitudinal 
polarization p o, i.e., the probability that the -D s i(2536) + helicity is zero. The other two 
probabilities are both equal to pn = p-1-1 = (1 — poo)/2. In the helicity formalism, the 
angular distribution in the decay chain _D sl (2536) + — > D* + K$, D* + — > D n + is given by 



d 3 N 
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cttcosa) d(3 <i(cos7) 4ir(l + 2Ra) 



x cos 7 



poo cos a + 



Poo . 2 
- — sin a 



+R A sin 2 7 



1 - Poo . 



sin (3 + cos /9(poo sin a + 



1 - Poo 



cos 2 a) 



+ 



^Ra(1-3poo) - 



sin 2a sin 27 cos {3 cos £ 



(3) 



The formula depends on three variables: p o, -Ra and £. Here V i?A = Aifi/A = z, 
where Ai >0 and A 0:0 are the helicity amplitudes corresponding to the D* + helicities ±1 
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and 0, respectively. They are related to S- and D-wave amplitudes in D s i(2536) + decay by 
= Ta^^T!^)' = \73^~ y/2D). Equation (j3J) allows one to extract poo and z from 
the D s i (2536) + angular distributions and to obtain D/S = y/2 (z— l)/(l+2z) = ^JTdJTs e ir > , 
where are the partial widths of -D s i(2536) + and 77 is the phase between D- and S- 
amplitudes. 

The last interference term in Eq. ([31 , with phase £, vanishes after integration over any 
angle. In particular, it does not appear in Ref. [§[ in the formulas for the two- and one- 
dimensional distributions d 2 N / d(pos a) oncost) and dN /cKposj). Therefore, in Refs. @, 0], 
only y/Rh = \z\ is measured for the _Di(2420) meson. This only constrains the possible 
ranges of Td/Ts and the phase 77. To determine them unambiguously, one needs to measure 
the phase £ and to fit the whole three-dimensional d 3 N/d((:os(y) dj3 d(cos7) distribution. 

The probability density function (PDF) for the unbinned maximum likelihood fit has the 
form 

V(a, (3, 7) = (1 - ft) ■ Tf r " C( "f' 7) + fb ■ Ptcfafr-Y). (4) 

a(cosa) dp d\posj) (e) avr 

It includes the efficiency corrections and the contribution of the background. The back- 
ground fraction /& = 528/6169 is estimated as the ratio of the number of entries in the 
sidebands and in the signal region, respectively. The signal and the sideband regions, de- 
fined as \AM D + I < 7 MeV/c 2 and 10 MeV/c 2 < \AM D + | <17 MeV/c 2 , respectively, are 

si si 

wider than in the D + tt~K + case since the background is lower. The PDF Vbck{&, P, 7), which 
is normalized to unity, is modelled using the sideband event distribution and the procedure 
described below. e(a,/3, 7) is the MC-determined efficiency The average efficiency (e) avr 
normalizes to unity the signal part of the PDF. It is recalculated in every iteration of the fit 
procedure as (e) avr = £i • I { « /// d(cosa f d p d{cosi f (a, P, 7) d(cosa) dp d(cosi). The sum is 
taken over 10x10x10 "bins" in a three-dimensional (cos a) x P x (cos 7) space. The efficiency 
map €i is determined from MC simulation, while the integral Ii of d 3 N/d(cosa) dp d(cos7) 
over each bin volume is calculated analytically. 

The density Vb c k(®, P, 7) of sideband events in the vicinity of (a,P,j) is calculated as 
follows. First, the three-dimensional (cos a) x/3x (cost) space is rescaled along each axis to 
the unit cube. This ensures that in the case of uniform distributions all three variables have 
the same "weight". Then, for the given point (a, P, 7), we find the volume V10 (Vn) of the 
smallest cube centered at this point and containing 10 (11) sideband events. The 10th (11th) 
event bisects a face of V w (Vn), so the half- weight of event 10 outside V w and the half- weight 
of event 11 inside Vn occupy the volume Vn — Vo- Therefore, we assign 10 full events to 
the volume V 10 , 5 = |(Vi + Vn) and estimate Vbck{a, P, 7) = 10/Vio.5/(47r ■ 528). Here, Air 
is the original volume of the (cosa)x/9x(cos7) space and 528 is the total number of sideband 
events. The resulting Vbck is thus normalized. To determine a systematic uncertainty due to 
this procedure, we use 20 or 50 points instead of 10. The changes are found to be negligible 
compared to statistical errors (see below). 

The advantage of this procedure is that for any one signal event Vbck is always determined 
from 10 (or 20, or 50) sideband events. Therefore the sideband fluctuations are much smaller 
than fluctuations of one signal event. This ensures the necessary degree of "smoothness" of 
the Vbck distribution. On the other hand, the typical volume V10.5 is about 10/528 ~ 0.02 of 
the whole (cos a) x/3x (cos 7) space volume. Vbck thus reproduces the background behaviour 
at this level of granularity. 

A similar method is used to construct the efficiency function e(a, P, 7). Due to the much 
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larger MC sample, instead of V10.5 we use the volume with 100 MC reconstructed events 
Vioo.5- We then determine the number of MC events generated there, N^, and calculate 
the efficiency as e(a, (3, 7) = lOO/iV™. As in the previous case, usage of the 50 or 200 
closest events instead of 100 reproduces the same results within the statistical errors, and is 
used to determine the systematic uncertainty of this method due to the efficiency 
The three-dimensional fit of all -D s i(2536) + signal entries to V(a,P,j) gives 

z = Aifl/Aofl = \Jr~kS = V3.6i0.3i0.lexp (±i • (1.27 i 0.15 i 0.05)). (5) 

Note that the angular distributions are sensitive only to cos£, not to £ itself. Therefore £ 
has a±( + 2nn ambiguity, and Ai }0 /A 0}0 is determined up to complex conjugation. The 
average _D s i(2536) + longitudinal polarization in the region xp>0.8 is measured to be poo — 
0.490 i 0.012 i 0.004. 

Systematic uncertainties are calculated as a sum in quadrature of the contributions listed 
in Table HTH MC simulation shows that the detector resolution in a, (3, 7 not only increases 
the final errors but also effectively decreases the parameter R\ by 0.13. The corresponding 
correction has already been applied to the above result. The systematic uncertainties in 
modelling Vbck and e are estimated by varying the number of closest points as explained 
above, and by using different sidebands: 7 MeV/c 2 < |AMd s1 (2536)+| < 10 MeV/c 2 plus 
17 MeV/c 2 < |AAf Dsl(2536) +| < 21 MeV/c 2 (instead of 10 MeV/c 2 < | AM Dal(2 53 6) + 1 < 
17 MeV/c 2 ). The errors due to statistical fluctuation of the MC sample, which is 15.4 times 

larger than data and has almost the same R\, £ and poo, are calculated as ^/l/15.4 = 0.25 
of the statistical errors. 



TABLE III: Systematic uncertainties for i?A, £ and pqq. 



Source 


Ra 


£ 


Poo 


Angular resolution 


0.05 


0.02 


0.001 


Modeling of Vbck and efficiency 


0.04 


0.00 


0.001 


Different sidebands 


0.03 


0.02 


0.002 


MC statistics 


0.07 


0.04 


0.003 


Total 


0.10 


0.05 


0.004 



The ratio of D- and S-wave amplitudes is found to be D / S = (0.63 i 0.07 i 0.02) • 
exp (±i • (0.76 i 0.03 i 0.01)). The relative phase is close to vr/4, (43.8 i 1.7 i 0.6)°. One 
can see that, contrary to the HQET prediction, the S-wave dominates. Its contribution to 
the total width is 1/(1 + \D/S\ 2 ) = T s /T total = 0.72 i 0.05 i 0.01. 

The background-subtracted efficiency corrected and normalized one-dimensional projec- 
tions of d 3 N / dfcos a) d(3 d(posj) distribution, together with the fit results, are shown in 
Fig. [5J The % 2 difference between the points and the projected fit results corresponds to 
a goodness-of-fit probability of about 60%. As mentioned earlier, one-dimensional projec- 
tions are not sensitive to the phase £. They are described instead by the following formulas 
derived from Eq. (j3J): 

" A 3 -{[(li J RA)i( J RA-l)p o]icos 2 a( J R A -l)(l-3poo)}, (6) 



dcosa 4(1 + 2/2/ 
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% = wTT2^o {[1 + 3jRa(1 " Poo)] + 2jRa(3/?00 _ x) cos2 (7) 

diV 3 



cicos7 2(l + 2i2/ 



■{cos 2 7 + J R A sin 2 7}. (8) 



0.08 
0.06 
0.04 
0.02 



"10 1 -10 1 

cosoc P cosy 

FIG. 5: Background subtracted, efficiency corrected and normalized one-dimensional cos a, (3, cos 7 
angular distributions. Projected results of the three-dimensional fit are shown by the solid curves. 

In spite of the complexity of Eq. (j3J), it depends only quadratically on V^a an d only 
linearly on cos£ and poo- The efficiency entering the PDF in Eq. (j3j) is almost constant in 
all three projections. When it is set to a constant value in the fit, the results change by less 
than 1/3 of the statistical error. The background fraction /& = 9% is small. Therefore one 
does not expect any significant biases of the fit results. 

To quantify this statement, 1000 samples of events are generated according to PDF 
Eq. (pi}. Each sample contains the same number of events as observed in data. The param- 
eters R\, £ and poo are set to the values determined from data. A three-dimensional fit is 
performed for each sample; it is verified that the fit results are not systematically biased and 
the errors are estimated correctly. The value of the overall likelihood function is measured 
to be worse than the one observed in data in 33% of cases. 

As a final check, the fit to data is repeated in different bins of the mass recoiling against 
the _D s i(2536) + , defined as ./(2E^ cam — E* + ) 2 — (p* + ) 2 , where all quantities are measured 

in the e + e~ center-of-mass frame. The parameters R\ and £ are found to be independent 
of -D s i(2536) + momentum or recoil mass within statistical errors. The recoil mass spectrum 
is shown in the top half of Fig. El The resolution is about 70 MeV/c 2 at 2 GeV/c 2 and is 
approximately inversely proportional to the recoil mass. There is an indication of two-body 
contributions from e + e" — > -D s i(2536) + A where X = D+, D* + and higher D** + resonances. 
This agrees with the _D s i(2536) + polarization spectrum, shown in the bottom half of Fig. [6l 
which also exhibits some structure at low recoil masses. This spectrum is obtained when R\ 
and £ are fixed to their values determined from the overall fit. At mass, as one expects, 
the longitudinal polarization is low. It then rises rapidly at the D* + mass and eventually 
reaches a plateau at poo ~ 0.5. 
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FIG. 6: 7Jgi(2536) + recoil mass spectrum (top), probability poo that L> s i(2536) + helicity is zero 
(bottom). The decrease at low values in the bottom plot can be attributed to the contribution of 
two-body e + e~ — > Z) s i(2536) + X reactions. 



CONCLUSIONS 

In conclusion, a new decay channel _D s i(2536) + — ► D + tt~K + is observed. The D + n~ pair 
is the only Dtt combination that cannot come from a D* resonance. It can be produced only 
through the virtual L>*°, broad D* ° or £^(2460)° resonances. In addition, the D + tt K + final 
state can be formed by two-body decays to a D + and a virtual K*° or higher K* resonance. 
No clear resonant substructure is found in the D + 7t~K + system. The ratio of branching 
fractions B(D+ X -> D+tx- K + ) f&(D+ x -> D*+K°) is measured to be (3.27 ± 0.18 ± 0.37)%. 

An angular analysis of the decay .D s i(2536) + — > D* + K$ is also performed. Since the c- 
quark is not infinitely heavy, HQET is violated and the _D s i(2536) + can contain an admixture 
of another J p = 1 + state with j = 1/2 and can decay in an S-wave. The energy release in 
this reaction is small. Therefore the D-wave is suppressed by the barrier factor {q/qo) 5 , 
where q is the relative momentum of D s i(2536) + decay products in the -D s i(2536) + rest 
frame, and go is a momentum scale characteristic of the decay. The S-wave contribution 
to the total width is proportional to q/qo and can be sizeable even if the mixing is small. 
Using an unbinned maximum likelihood fit to the three angles in the D sl (253Q) + -> D* + K° s , 
D* + — > D°7i + decay chain, we measure (up to a complex conjugation) the ratio of S- and 
D-wave amplitudes: D/S = (0.63±0.07±0.02) -exp (±i ■ (0.76 ± 0.03 ± 0.01)). The S-wave 
dominates, and its contribution to the total width is Ts/T tota i = 0.72 ± 0.05 ± 0.01. This 
result allows to calculate the mixing angle in the theoretical models with a known value of 
parameter q 0, EH . 

The spin of high momentum D s i(253Q) + mesons produced in e + e~ annihilation prefers 
to align transversely to the momentum. The probability that a _D s i(2536) + with xp> 0.8 
has zero helicity is found to be poo = 0.490 ± 0.012 ± 0.004. Assuming the HQET relation 
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Poo = |(1 — W3/2)) \1M this implies a value of the Falk-Peskin parameter, w 3 / 2 = 0.266 ± 
0.018 ± 0.006, in this momentum region. This value is close to the prediction of Ref. 16], 
W3/2 ~ 0.254, obtained for the entire momentum region, although the applicability of the 
perturbative QCD fragmentation model for D** + mesons is questionable. 
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